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Predictions on how animals respond to habitat changes 
are the primary aim of many conservation studies. De¬ 
velopment of easy wildlife habitat models is an important 
tool for conservation and ecosystem management (Gon¬ 
zalez et al. 1992, Donazar et al. 1993). Progress has been 
made using Generalized Linear Models (GLMs) (Dobson 
1983, McCullagh and Nelder 1989) to summarize the re¬ 
lationships between species distributions and environ¬ 
mental variables (Vincent and Haworth 1983, Nicholls 
1989, Donazar et al. 1993). 

It is known that patterns and processes in nature are 
sensitive to the scale at which they are viewed (Cody 


1985, Wiens et al. 1987, Wiens 1989, Levin 1992, Lima 
and Zollner 1996). The scale at which systems are studied 
has a powerful influence on final conclusions and spe¬ 
cies-habitat relationships determined at one scale may 
not apply to others. Populations are influenced by the 
complex arrangement of habitat patches within land¬ 
scapes and multiscaled studies seem to be the proper way 
to approach their study (Wiens 1989, Levin 1992). 

The Golden Eagle (Aquila chrysaetos) is a raptor with a 
widespread distribution in the northern hemisphere. In 
North America, Steenhof et al. (1997) showed an impor¬ 
tant interaction between jackrabbit (Lepus califomicus) 
abundance and weather on eagle reproduction and more 
recent work using radiotracking data (Marzluff et al. 
1997) has noted the preference of Golden Eagles for 
some habitat types, particularly shrub and open lands. In 
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Europe, different qualitative descriptions around nest 
sites have been published (Tjernberg 1983, Watson 
1997), and McGrady et al. (1997) constructed a model 
that delineated the area over which eagle pairs range and 
habitats of particular importance. More general ap¬ 
proaches should be used as comparative parameters to 
consider entire populations and to reduce individual var¬ 
iability (White and Garrot 1990, Aebischer et al. 1993). 

The purpose of this paper is to make a mathematical 
description of Golden Eagle breeding sites in southeast¬ 
ern Spain, where one of the highest densities of this spe¬ 
cies has been reported (Sanchez-Zapata et al. 1995). We 
propose a two-scale approach, considering at first re¬ 
sponses around nest sites, and then a landscape-level 
analysis to evaluate the influence of the matrix on breed¬ 
ing territories. 

Study Area and Methods 

The study area covered the Murcia region, a 11317 
km 2 area located in southeastern Spain with numerous 
mountains ranging from 0-2000 m elevation. The cli¬ 
mate is Mediterranean arid and semiarid with a mean 
annual rainfall of 300 mm. Vegetation has a mosaic struc¬ 
ture with cultivated lands (54%), grasslands and shrub- 
lands (28%), forest (15%) and open lands (3%) (Alcaraz 
et al. 1991). 

All the territories known to be occupied by Golden 
Eagles at least once during the period 1985-97 were con¬ 
sidered (Sanchez-Zapata et al. 1995). The location of 
breeding territories was incorporated into a Geographic 
Information System (IDRISI, Eastman 1992) using the 
UTM grid of 1 km 2 cells. For the first small-scale land¬ 
scape approach the 1 km 2 cells were aggregated into 9 
km 2 (3 X 3 km) cells, so the regional map of 11 317 km 2 
cells was transformed into a map with 1381 cells of 9 km 2 . 
The large-scale landscape analysis was focused on 88 cells 
of 100 km 2 (10 X 10 km). 

The same GIS was used to characterize the breeding 
sites using the following variables (Table 1): (1) SLOPE 
and LAND USE—slope (° from horizontal) was calculat¬ 
ed from a Digitized Land Model 1:100000 (Servicio Car- 
tografico Espahol) by comparing the altitude of each ba¬ 
sic cell (200 X 200 m) with that of neighboring cells to 
the north, south, east and west. An average value for the 
different 200 X 200 m cells was calculated. Slope for larg¬ 
er cells (3X3 km and 10 X 10 km) was obtained as the 
mean value of 200 X 200 m subcells. These values ranged 
from 0-24.2 at the 9 km 2 scale and from 0.2-13.5 at the 
100 km 2 scale. Different land-use classes were obtained 
from maps of the Ministerio de Agricultura (1:200 000) 
as proportions of cell area (9 km 2 and 100 km 2 ) covered 
by each. New categories were formed by combining re¬ 
lated land-use cover categories (e.g., lemon, orange and 
other fruit trees were combined to give a single arbore¬ 
ous intensive agriculture category). (2) EDGE—edge was 
measured as the length (km) of edges between different 
land uses using the digitalized land-use map and ATLAS 
GIS software. (3) STRUCTURE—number and size (ha) 
of the different patches of natural vegetation obtained 
from maps of the Direccion General de Produccion 
Agraria (1:200000). 


Table 1. Variables used to characterize the breeding ar¬ 
eas of Golden Eagles in southeastern Spain. 


Variables Used in General Linear Model 
Land Use Categories 

AINTA—% of cell covered by arboreous intensive ag¬ 
riculture, such as lemon and orange trees. 

HINTA—% of cell covered by herbaceous intensive ag¬ 
riculture, such as vegetable crops. 

AEXTA—% of cell covered by arboreous extensive ag¬ 
riculture, such as olive and almond trees. 

HEXTA—% of cell covered by herbaceous extensive ag¬ 
riculture, such as cereal crops. 

SHRUB—% of cell covered by shrubland. 

FOREST—% of cell covered by forest, mainly Pinus ha- 
lepensis. 

SHF—% of cell covered by mixed shrubforest. 

SLOPE—topographic irregularity index. 

Edges 

EAEA—length (km) of edges between intensive and ex¬ 
tensive agriculture. 

FOIA—length (km) of edges between intensive agri¬ 
culture and forest. 

IASH-length (km) of edges between intensive agricul¬ 
ture and shrubland. 

IASF—length (km) of edges between intensive agricul¬ 
ture and mixed shrubforest. 

FOEA—length (km) of edges between forest and ex¬ 
tensive agriculture. 

EASH—length (km) of edges between extensive agri¬ 
culture and shrubland. 

EASF—length (km) of edges between extensive agri¬ 
culture and mixed shrubforest. 

FOSH—length (km) of edges between forest and 
shrubland. 

FOSF—length (km) of edges between forest and mixed 
shrubforest. 

SHSF—length (km) of edges between shrubland and 
mixed shrubforest. 

Structure 

PATCH—number of land-use patches per cell. 

RICHNESS—number of different land-use patches per 
cell. 

DIVERSITY—diversity (Shannon-Weiner) of land uses 

NFOREST—number of forest patches per cell. 

SFOREST—mean size (ha) of forest patches per cell 

NSHRUB—number of shrubland patches per cell. 

SSHRUB—mean size (ha) of shrubland patches per 
cell. 

NSHF—number of mixed shrub-forest patches per cell 

SSHF—mean size (ha) of mixed shrub-forest patches 
per cell. 

NNAT—number of natural vegetation patches per cell. 

SNAT—mean size (ha) of natural vegetation patches 
per cell. 
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Table 2. Response of Golden Eagles to the different habitat variables considered in southeastern Spain. % dev: 
deviance explained (ns—not significant, * P < 0.05, ** P < 0.01, *** P < 0.001. Responses: + s-shaped function, + + 
bell-shaped function, — s-shaped function. 




Scale 



3X3 


10 X 10 


% Dev 

Response 

% Dev 

Response 

Land Use 





AINTA 

ns 


ns 


HINTA 

9.18* 

— 

11.37* 

— 

AEXTA 

ns 


ns 


HEXTA 

7.65** 

— 

ns 


SHRUB 

ns 


ns 


FOREST 

13.09*** 

+ 

22-81*** 

+ 

SHF 

7.60* 

+ + 

21.88*** 

+ 

SLOPE 

28.89** 

+ + 

4212*** 

+ 

Edges 





EALA 

6.76** 

— 

9.61** 

— 

FOLA 

ns 


ns 


IASH 

ns 


ns 


IASF 

ns 


ns 


FOEA 

5.46* 

+ 

14.45* 

+ + 

EASH 

3.98** 

— 

ns 


EASF 

ns 


ns 


FOSH 

2.88** 

+ 

ns 


FOSF 

4.95* 

+ + 

26.94*** 

+ 

SHSF 

9.59*** 

+ 

ns 


Structure 





PATCH 

2.58** 

— 

ns 


RICHNESS 

3.11* 

+ + 

ns 


DIVERSITY 

1.88* 

+ + 

7.33* 

+ + 

NFOREST 

7 80*** 

+ + 

ns 


SFOREST 

11.34*** 

+ + 

17.84* 

+ + 

NSHRUB 

ns 


ns 


SSHRUB 

ns 


ns 


NSHF 

9.80** 

+ + 

8.31** 

+ 

SSHF 

5.24* 

+ + 

ns 


NNAT 

1.96** 

+ 

ns 


SNAT 

11.60*** 

+ + 

15.64* 

+ + 


Due to the effect of increasing sizes of cells, many 
of them included large areas of sea and adjacent re¬ 
gions that were not censused. Therefore, these cells 
were excluded from the data analysis making fewer 
breeding territories in the large-scale study (76 vs. 40 
territories). 

We used Generalized Linear Models (GLMs) to con¬ 
struct models of the breeding density of Golden Eagles 
(Dobson 1983, McCullagh and Nelder 1989, Nicholls 
1989). For density response variables (number of breed¬ 
ing territories), the Poisson distribution was an adequate 
error function (Vincent and Haworth 1983) and the dis¬ 
crete Poisson function an appropiate link function (L = 
e (a+bi Xl +.. ,+b k x k )y This meant that the number of breeding 


sites in an area was a discrete, s-shaped function when 
the linear predictor was the first order polynomial or a 
bell-shaped function for second order polynomials (San- 
chez-Zapata and Calvo 1999). 

For regression analysis, we used the program STATIS- 
TIX (Analytical Software 1992) following a forward step¬ 
wise analysis (Donazar et al. 1993). Each explanatory var¬ 
iable was tested for significance in turn. The variable 
contributing to the largest significant change in deviance 
from the null model was then selected and fit to the mod¬ 
el. Once a variable was fit to the model, we tested if the 
addition of a second variable significantly improved the 
model. We chose a 5% level of significance to include a 
variable in a model. 
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Table 3. General Linear Model for Golden Eagle breeding density in southeastern Spain, % dev: deviance explained. 



Coefficient 

SE 

P 

% Dev 

9 km 2 Scale 

Constant 

-6.55778 

0.62282 

0.0000 

30.16 

SLOPE 

0.60749 

0.12321 

0.0000 


SLOPE 2 

-0.0167 

0.00562 

0.0000 


SHSF 

2.034e~ 4 

8.296e -5 

0.0142 


Constant 

-6.64018 

0.63324 

0.0000 

30.32 

SLOPE 

0.60065 

0.12638 

0.0000 


SLOPE 2 

-0.01693 

0.0058 

0.0035 


SHF 

0.23580 

0.10192 

0.0207 


Constant 

-6.86474 

0.65345 

0.0000 

31.43 

SLOPE 

0.54579 

0.13264 

0.0000 


SLOPE 2 

-0.01445 

0.00601 

0.0162 


SSHF 

8.415e 7 

3.826e- 7 

0.0279 


SSHF 2 

-1.447e- 13 

7.055e -14 

0.0402 


100 km 2 Scale 

Constant 

-0.75439 

0.72879 

0.0000 

43.70 

SLOPE 

0.38313 

0.05857 

0.0000 


SHRUB 

0.06663 

0.03237 

0.0396 



Results and Discussion 

SLOPE was the most important variable at both scales 
explaining a higher percentage at the 100 km 2 scale 
(28.89% and 42-12%, respectively). When scale changed, 
the response of eagles changed from quadratic to linear. 
Pinus halepensis forests (FOREST) was the second most 
important variable (13.09% and 22.81% for each scale) 
with a similar percentage of deviance explained by mixed 
shrubforest at the 10 X 10 km scale (21.88%). Intensive 
agriculture was negatively associated with eagles (HINTA 
9.18 and 11.37%, respectively), while cereal crops were 
negatively correlated only at the smaller scale (HEXTA 
7.65%) (Table 2). 

Edges between land uses were also important at the 
larger scale, but explained low percentages of deviance 
at the 3X3 km scale. Edges between forest and mixed 
shrubforest (FOSF 26.94%) and edges between forest 
and extensive agriculture (FOEA 14.45%) were the most 
explanatory edge variable. Eagles responded negatively 
to edges between extensive and intensive agriculture at 
both scales (EAIA 6.76% and 9.61%, respectively). There 
was a negative relationship between eagle densities and 
edge between shrub and extensive agriculture only at the 
9 km 2 scale (EASH 3.98%) (Table 2). 

Percentages of deviance explained by landscape struc¬ 
ture were generally low, except for natural vegetation and 
forest patch sizes (SNAT 11.60%, 15.64% and SFOREST 
11.34%, 17.84% for each scale). At the 9 km 2 scale, the 
number of patches of natural vegetation and eagle den¬ 
sities seemed to be negatively related (PATCH 2.58%) 
(Table 2). 

Because slope accounted for the higher percentages of 


deviance, models were constructed entering SLOPE as 
the first variable. Only a small reduction in deviance was 
obtained by including other variables. At the larger scale, 
the model was more explanatory than at the 3X3 km 
scale (Table 3), 

Because most Golden Eagles in Murcia nest on cliffs, 
slope was the most important variable in predicting its 
breeding densities. The linear response at the larger 
scale suggested that eagles preferred the bigger moun¬ 
tain systems of the region. 

The primary factor influencing Golden Eagle breeding 
success is food availability (Steenhof et al. 1997). Eagles 
prey on medium-sized mammals such as jackrabbits (Le- 
pus spp.) and rabbits (Oryctolagus cuniculus ) (Steenhof et 
al. 1988, 1997, Watson 1997), which are very common in 
shrublands in Mediterranean areas (Moreno and Villa- 
fuerte 1995, Palomares and Delibes 1997). These open 
lands, where vegetation structure favors prey detection 
and hunting success (Tjernberg 1983, Marzluff et al. 
1997, McGrady et al. 1997), were the second most im¬ 
portant variable after slope in modeling breeding densi¬ 
ties. Although forests seemed to be an important factor, 
this could have been a consequence of their distribution 
in the mountain systems of the area (Chaparro 1996). 

Eagle densities were negatively correlated with irrigat¬ 
ed crops, possibly because of the high number of people 
working in fields and their low prey populations. The 
increase in power lines associated with irrigated land 
could have also been an important negative factor for 
Golden Eagles, as electrocution is the main cause of mor¬ 
tality for many eagle species (Gonzalez et al. 1990, Ferrer 
and Hiraldo 1992, Sanchez-Zapata et al. 1995). 
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Resumen. —Utilizando Generalized Linear Models 
(GLMs) examinamos la densidad reproductiva del Aguila 
real Aquila chrysaetos en relation con los usos del suelo, 
los hordes entre usos y la estructura del paisaje en el 
sureste de Espana. Las respuestas se compararon a dos 
escalas de paisaje. La pendiente fue la variable mas im- 
portante para predecir la densidad reproductiva. Las 
manchas de vegetation natural se relacionaron positiva- 
mente con las aguilas mientras que la agricultura inten- 
siva se correlation© de manera negativa. El matorral pa- 
rece tener efectos positivos importantes, probablemente 
al incrementar la disponibilidad de alimento. 

[Traduction de Autores] 

Acknowledgments 

We would like to thank Sergio Eguia, Jose E. Martinez, 
Miguel A. Sanchez and Javier Royo for their assistance in 
the fieldwork. 

Literature Cited 

Aebischer, N.J., RA. Robertson and RE. Kenward. 
1993. Compositional analysis of habitat use from ani¬ 
mal radio-tracking data. Ecology 74:1313-1325. 
Alcaraz, F., P. Sanchez-Gomez, A. de la Torre, S. Rios 
and J. Alvarez Rogel. 1991. Datos sobre la vegeta¬ 
cion de Murcia (Espana). Guia Geobotanica de la Ex¬ 
cursion de las XI Jornadas de Fitosociologia. DM- 
PPU. Madrid, Spain. 

Analytical Software. 1992. STATISTIX. St. Paul, MN 
U.S.A. 

Chaparro, J. 1996. Distribution potencial del bosque y 
de sus especies arboreas en zonas mediterraneas semi- 
aridas: modelos y aplicaciones. Ph.D. dissertation, 
Universidad de Murcia, Murcia, Spain. 

Cody, M.L. [Ed.]. 1985. Habitat selection in birds. Aca¬ 
demic Press, San Diego, CA U.S.A. 

Dobson, A.J. 1983. Introduction to statistical modeling. 

Chapman and Hall, London, U.K. 

Donazar, J.A., F. Hiraldo and J. Bustamante. 1993. Fac¬ 
tors influencing nest site selection, breeding density 
and breeding success in the Bearded Vulture ( Gypae- 
tus barbatus). J. Appl. Ecol. 30:504-514. 

Eastman, J.R. 1992. IDRISI 4.0. Worcester, MA U.S.A. 
Ferrer, M. and F. Hiraldo. 1992. Man-induced sex-bi¬ 
ased mortality in the Spanish Imperial Eagle. Biol. 
Conserv. 60:57-60. 

Gonzalez, L.M., J. Bustamante and F. Hiraldo. 1990. 
Factors influencing the present distribution of the 
Spanish Imperial Eagle Aquila adalberti. Biol. Conserv. 
51:311-319. 

-,-and -. 1992. Nesting habitat selec¬ 
tion by the Spanish Imperial eagle Aquila adalberti. 
Biol. Conserv. 59:45-50. 


Levin, S.A. 1992. The problem of pattern and scale in 
ecology. Ecology 73:1943-1967. 

Lima, S.L. and P.A. Zollner. 1996. Towards a behavioral 
ecology of ecological landscapes. TREE 11:131-135. 

Marzluff, J.M., S.T. Knick, M.S. Vekasy, L.S. Schueck 
AND T.J. Zarriello. 1997. Spatial use and habitat se¬ 
lection of Golden Eagles in southwestern Idaho. Auk 
114:673-686. 

McCullagh, P. and J.A. Nelder. 1989. Generalized lin¬ 
ear modeling. Chapman and Hall, London, U.K. 

McGrady, M.J., D.R.A. McLeod, S.J. Petty, J.R. Grant 
and I.P Bainbridge. 1997. Golden Eagles and forest¬ 
ry, Res. Inform. Note 292. Forestry Commission Re¬ 
search Agency. Aberdeen, U.K. 

Moreno, S. AND R. Villafuerte. 1995. Traditional man¬ 
agement of shrubland for the conservation of rabbits 
Oryctolagus cuniculus and their predators in Donana 
National Park, Spain. Biol. Conserv. 73:81-85. 

NlCHOLLS, A.O. 1989. How to make biological survey go 
further with generalized linear models. Biol. Conserv. 
50:51-75. 

Palomares, F. and M. Delibes. 1997. Predation upon Eu¬ 
ropean rabbits and their use of open and closed 
patches in Mediterranean habitats. Oihos 80:407-410. 

Sanchez-Zapata, J.A., M.A. Sanchez, J.F. Calvo and M.A. 
Esteve. 1995. Ecologia de las aves de presa en la re¬ 
gion de Murcia. Cuadernos de Ecologia y Medioam- 
biente. Universidad de Murcia, Murcia, Spain. 

-and J.F. Calvo. 1999. Raptor distribution in rela¬ 
tion to landscape composition in semi-arid Mediter¬ 
ranean habitats./. Appl. Ecol. 36:254-262. 

Steenhof, K., M.N. Kochert and J.H. Doremus. 1988. 
Nesting of subadult Golden Eagles in southwestern 
Idaho. Auk 100:743-747. 

-,-AND T.L. McDonald. 1997. Interactive ef¬ 
fects of prey and weather on Golden Eagle reproduc¬ 
tion. /. Anim. Ecol. 66:350-362. 

TjERNBERG, M. 1983. Population density of Golden Eagle 
in relation to nest-size and food availability. In Breed¬ 
ing ecology of the Golden Eagle, Aquila chrysaetos (L,), 
in Sweden. Report 10. Swedish Univ. Agric. Sci., Dept. 
Wildlife Ecology. Uppsala, Sweden. 

Vincent, J. and M. Haworth. 1983. Poisson regression 
models of species abundance./. Biogeogr. 10:153-160. 

Watson, J. 1997. The Golden Eagle. T. & A.D. Poyser, 
London, U.K. 

White, G.C. and R.A. Garrot. 1990. Analysis of wildlife 
radio-tracking data. Academic Press, London, U.K. 

Wiens, J.A. 1989. Spatial scaling in ecology. Fund. Ecol. 3: 
385-397. 

-, J.T. Rotenberry and B. van Horne. 1987. Habi¬ 
tat occupancy patterns of North American shrubsteppe 
birds: the effects of spatial scale. Oikos 48:132-147. 

Received 25 February 1999; accepted 28 October 1999 









